A dynamic scheme for generating number squeezing in Bose-Einstein condensates 

through nonhnear interactions 
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We develop a scheme to generate number squeezing in a Bose-Einstein condensate by utilizing 
interference between two hyperfine levels and nonlinear atomic interactions. We describe the scheme 
using a multimode quantum field model and find agreement with a simple analytic model in certain 
regimes. We demonstrate that the scheme gives strong squeezing for realistic choices of parameters 
and atomic species. The number squeezing can result in noise well below the quantum limit, even 
if the initial noise on the system is classical and much greater than that of a poisson distribution. 



Introduction. — The experimental realization of Bose- 
Einstein condensates (BECs) has allowed the creation 
of macroscopic quantum systems that are highly con- 
trollable, and hence provide an excellent system to test 
predictions of many body quantum dynamics. The gen- 
eration of nonclassical states in BECs, such as number 
squeezed states, would allow for measurements of parti- 
cle number statistics that differ from classical predictions 
However, while the motional state of the atoms is 
reasonably simple to manipulate, the quantum statistics 
referring to the number of particles present is difficult 
to control, with BECs typically produced with a large 
(5%) shot to shot variation in the number. The genera- 
tion of nonclassical states in BECs is currently of great 
interest [H, as it could potentially enhance the sensitiv- 
ity of atomic interferometers [3[ used to measure electric, 
magnetic, gravitational fields, accelerations, and atomic 
interactions. The generation of nonclassical states also 
provides a method for testingthe fidelity of recent quan- 
tum state transfer schemes [j], and an atom laser pro- 
duced from a number squeezed BEC will have a reduced 
linewidth [5]. 

The generation of nonclassical states via self- 
interaction in samples of cold atoms has been considered 
before 0, 0, i, 0, M [Hi. The schemes proposed in [6] 
and describe the generation of quadrature squeezing, 
which requires the use of a well-defined phase reference 
in order to be observed. A well-defined phase reference 
is difficult to obtain in atom optical systems, especially 
in the presence of strong nonlinear it ies, which are re- 
quired to produce the quadrature squeezing, thus mak- 
ing the schemes somewhat unrealistic. In addition, both 
of these schemes assume the BEC is initially in a coher- 
ent state, rather than a more realistic statistical mixture 
of coherent states with random phases. 0, 0, [l^ have 
demonstrated that the atomic nonlinearity can be used 
to generate number difference squeezing, creating a state 
with angular momentum projection below the standard 
quantum limit (spin squeezing). These schemes assume 
that the total number of particles is initially well defined. 
Recently, Esteve et al. [2] have directly observed squeez- 



ing in the number difference between two adjacent lattice 
sites. 

In this Letter we describe a scheme that allows the 
creation of absolute number squeezing in a BEC (as 
opposed to number difference squeezing), which is ex- 
perimentally realistic, utilizes only the relatively simple 
experimental technique of Ramsey interferometery, and 
does not require manipulation of the scattering length 
via a Feshbach resonance or a coherent phase reference 
for the atoms. We show that our scheme achieves num- 
ber squeezing below the quantum limit even if there is 
initially considerable classical noise on the system, and 
that this result holds even when we realistically assume 
that the initial state of the BEC is a statistical mixture 
of states with random phase. 

We consider a BEC with two internal states confined 
to an optical trap, with all the atoms initially in one 
state. A short state-dependent coupling is applied, trans- 
ferring a small fraction of the population to another state. 
The system is then left to evolve for some time, allow- 
ing nonlinear interactions and interference between the 
two states, before the coupling is applied for a second 
time, transferring some of the population back to the 
initial state. Provided the 5- wave scattering lengths of 
the atoms in the different internal states are not all iden- 
tical, and by choosing appropriate coupling strengths, 
hold times and trap geometry, it is possible to gener- 
ate number squeezing. An important difference between 
our scheme and the one demonstrated by Esteve et al. 
[2] is that our scheme is based on dynamic interfer- 
ence between the two modes to obtain absolute num- 
ber squeezing in one of the modes, whereas Esteve et 
al. obtain their squeezing by adiabatically changing the 
potential to one where the ground state of the system 
exhibits number difference squeezing. Scheme. — Our 
proposed scheme is based on a Ramsey interference ex- 
periment between two hyperfine states of sodium, namely 
\F = l,mF = +1) = |1) and \F = 2,mF = 0) = |2). The 
timing for our scheme is outlined in Figure [H We begin 
with all the atoms in a BEC in the |1) state in an optical 
trap. At t = to, the microwave coupling is turned on for 
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FIG. 1: Timing for the coupling pulses in the proposed ex- 
periment. The coupling field is turned on at to, and then 
off again at ti. After a duration thoid, the coupling pulse is 
turned back on again at t2 and finally turned off at ta. After 
this, the population of state |2) atoms is measured. 

a brief duration of time Ati = ti — to. During this time, 
a fraction of the atoms are then transferred to state |2). 
The coupling is switched off, and the system is left to 
evolve for an amount of time thoid = ^2 — ^i, before we 
interfere the two modes with a second microwave cou- 
pling pulse for a duration At2 = ^3 — ^2- We assume that 
this pulse is phase locked to the first, with an adjustable 
phase shift 0. This phase shift is an important ingredi- 
ent, as it is what allows us to ensure that the two modes 
interfere in such a way as to produce number squeezing, 
independently of thoid, which controls the depth of squeez- 
ing. Nandi et al [11] have recently published a scheme 
based on a Ramsey-Borde interferometer, but without 
including this adjustable phase shift (j) at the second cou- 
pling pulse. As a result, their squeezing is only observed 
at particular values of thoid, when the nonlinear phase 
shift acquired during the hold time is appropriate. The 
inclusion of this phase shift gives us an additional degree 
of freedom in our system, as we can choose thoid inde- 
pendently of the phase shift required to produce number 
squeezing, and thus optimize our squeezing depth. Af- 
ter ^3, we separate the two modes with a magnetic field, 
and count the number of atoms in mode |2) to determine 
the number statistics. The Hamiltonian for the system 
is = Ho^Hcit), with 

+ E ^ /v^Iw^jwv^.wv^.w (1) 

and 

= J (fil^(t)e^^V^J(r)^i(r)+h.c.) d^r, (2) 

where '0i(r) represents the annihilation operator for state 
Hj = ^^V^ + Kpt(r) + (j — 1)M is the single particle 
Hamiltonian, and lopt(r) is the optical dipole potential. 
Q{t) represents the microwave coupling field, which is 
switched on and off to control the coupling between the 
two hyperfine levels. The phase of this RF field can also 
be tuned between each pulse. We will assume that the 
coupling is on resonance, such that rt{t) = ^]oe~*'^^ with 
hd the hyperfine splitting between |1) and |2), and is 
the Rabi frequency. 



Analytic model — We first consider a two mode 
model, which demonstrates how atomic nonlinearities 
can be used to generate number squeezing. A two-mode 
model can be derived from Eq. ([1]) and Eq. ([2]) by assum- 
ing that the atoms remain in the ground motional state 
of the optical trap. With this assumption, the modified 
Hamiltonian for the system is = Ho + TYc, with 

Ho = hS 0^20.2 + E hxijCi\aid}-aj (3) 

He = h(n{t)e''f'alai^h.c}j (4) 

where ai {0.2) annihilates an atom from state |1) (|2)), 
and Xij = / IV^olr)!^ <i^r, where '0o(r) is the ground 
state wavefunction of the optical potential. 

Beginning with the coupling initially switched off 
{Qo = 0), and assuming that our initial state is a Poisso- 
nian mixture of number states for mode |1), (or, equiva- 
lently, a mixture of coherent states with random phase) 
and vacuum for mode |2) we notice that the evolution is 
trivial, as the density operator commutes with Ho. At 
to, the coupling is turned on for a duration Ati cou- 
pling a fraction of the atoms (sin^(r^oAti)) into mode 
1 2), creating a relative phase between the two modes. 
This relative phase is then 'sheared' during thoid by evo- 
lution under Ho due to the Kerr effect. We note that if 
Xii = X22 = X12, there is no phase shearing due to this 
effect, and our scheme does not work. We chose sodium 
as our atomic species, as it has a relatively large differ- 
ence between the scattering lengths of|F = l,mF = +l) 
and \F = 2,mF = 0). If we assume that Ati is suffi- 
ciently short that we can ignore the effect of Ho during 
the coupling time, the density matrix for the system after 

^hold is 

ni ,n2 ,mi=0 

X \ni,n2){mi,ni n2 - mi\ , (5) 

with <l>ni,n2,mi Xll(^l(^l " 1) " mi(mi - 

1)) + X22((^l - nl)(2n2 + ni - mi - 1)) + 
2xi2(nin2 - mi(ni + n2 - mi)) and An^^mi,n2 = 

^y ni!mi!n2!(ni+n2— mi)! 

a{ti) = aocos^i, /3{ti) = — mosin^i, where Oi = Qo^ti 
and |aoP = A/'o, the mean number of atoms initially 
in state |1). We note that at this point, both modes 
still contain a Poissonian number distribution. Fi- 
nally, we describe the dynamics caused by the second 
microwave pulse in the Heisenberg picture, by not- 
ing that the Heisenberg operators after the second 
pulse are clih = ai(0)cos^2 — ^2(0) sin^2e*'^, and 
= a2(0)cos6>2 — mi(0) sin 6^26"*"^, with O2 = l^oAt2, 
where At2 is the duration of the second microwave pulse, 
and (j) is the phase of the microwave field relative to the 
first pulse. Again, we have assumed that the duration 
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of the pulse is sufficiently short that we can ignore 
the evolution due to Hq. As we are only interested 
in the number statistics, we can neglect the rest of 
the evolution after the second microwave pulse, as the 
number operators for both modes commute with TYo- 
Assuming we can distinguish state |2) atoms from state 
|1) atoms, we define the normalized number variance 
for state |2) atoms as v{N2) = ((iV|) - (/^2)^)/(^2), 
N2 = CLljja2H' Figure [2] shows the v{N2) as a function 
of thoid and 0, using the scattering properties of sodium 
(ail = CLi2 = 2.8 nm, a22 = 3.0 nm [lH in a 500 Hz 
spherical harmonic trap. For some values of thoid and 0, 
v{N2) dips below 0.01, as compared to the quantum limit 
v{N2) = 1 associated with a coherent state, indicating 
significant number squeezing. The parameter space for 
this model is quite large, as we can adjust the length 
of the first and second coupling pulses, the hold time, 
and the phase of the second coupling pulse. If we did 
not have the ability to adjust the phase of the second 
coupling pulse, we would be constrained to a vertical 
line in figure [21 and not necessarily be able to access 
the optimum value of the squeezing. We found that 
the best number squeezing was obtained when the first 
coupling pulse was quite weak (approximately 8% of 
the atoms transferred). We also found we could still 
get a good level of squeezing when we began with an 
initial state which had number fiuctuations 150 times 
larger than a Poissonian distribution (about 5% shot 
to shot fiuctuations in the number) (see Figure [3]). 
When starting with a such an initial condition, the best 
squeezing was found when the first beam splitter was 
relatively weak. This is due to the fact that the addition 
of vacuum to a super-Poissonian number distribution 
drives it to a Poissonian number distribution. 

Multimode model. — To investigate if the approxima- 
tions we made in the previous section were valid, we per- 
formed a ID multimode simulation of the system using 
a stochastic phase space method. Specifically, we uti- 
lize a Truncated Wigner (TW) approach [13] • We reduce 
Eqs. ([1]) and ([2]) to one dimension by integrating out the 
dynamics in the y and z dimensions. A Fokker Plank 
equation (FPE) is then found from the master equation 
for the system using the Wigner representation. This 
equation can then be converted into a set of stochas- 
tic partial differential equations (SPDEs), which can be 
solved numerically. By averaging over many trajecto- 
ries with different noises, expectation values of quanti- 
ties corresponding to operators in the full quantum field 
theory can be extracted. When converting our FPE to 
a SPDEs, we ignore third and higher order derivatives 
in the FPE, as these terms do not have a simple map- 
ping to the stochastic PDEs, and can be assumed to be 
negligible when the field has a high occupation number 
[l3|. This truncated Wigner approximation will eventu- 
ally fail, as it can not describe negative components of 
the Wigner function, which eventually occur when evolv- 
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FIG. 2: (color online) \ogiQ(v{N2)) as a function of thoid and 
(/). For some values of thoid and 0, v{N2) dips below 0.01, 
indicating significant squeezing. Parameters: xii — X12 = 
0.018 s~\ X22 = 0.019 s~\ The strength of the coupling 
pulses was Oi = 0.3 rad and O2 = 0.025 rad for the first and 
second coupling pulses respectively. The initial occupation of 
the BEG was chosen to be a Poisson distribution with (N) = 
10^. These parameters correspond to the scattering properties 
of sodium in a 500 Hz spherical harmonic trap. 



ing under a Hamiltonian such as Eq. ([T|). However, we 
have checked our simulations in limits where the multi- 
mode dynamics can be neglected and they agree with the 
results obtained from our two-mode analytical model. In 
addition, over our simulation times no anomalous results 
such as signficant negative densities were seen, indicating 
that the truncation of third order derivatives was a valid 
approximation. The SPDEs describing the one dimen- 
sional system are 

ihipiix) = Ciipi{x) ^ hn{t)ip2{x) (6) 

ihMx) = C2Mx)^^^%t)Mx). (7) 



with 



+ Ujjili^jf - 1/dx) + Uijm^ - l/{2dx)) , (8) 



where dx is the grid spacing of the numerical simulations. 
The terms inversely proportional to dx compensate for 
the mean field of the vacuum, which is nonzero in the 
Wigner approach. The noise on the initial conditions for 
each trajectory of the evolution of these equations was 
chosen such that they corresponded to the specific initial 
state of interest. 

Figure [3] shows the results from the TW simulation 
for two different parameter regimes. The left column 
shows a case where there is excellent agreement between 
the multimode model and the two-mode analytic model. 
However, the right column shows a case where the com- 
parison between the two models is poor, due to signifi- 
cant multimode dynamics preventing the system acting 
as a two mode system. The multimode dynamics is a 
consequence of the unequal scattering lengths, meaning 
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FIG. 3: (color online) Results from the multimode TW model 
compared to the analytic two mode model, for two differ- 
ent parameter regimes (left column and right column respec- 
tively), (a) and (b) show the density profile of state |2) atoms 
at t2 (blue solid line), compared to the density profile at ti. 
(c) and (d) show number variance v{N2) for an initial state 
that is a mixture of coherent states (solid black line) and a 
mixture of states with 5% number uncertainty (dashed line), 
while (e) and (f) show (iV2), at ts. In (c) and (e), there is ex- 
cellent agreement between the multimode TW model (red cir- 
cles) and the two mode model (solid black trace). In (d) and 
(f), there is significant disagreement between the two models 
because a larger fraction of atoms is transferred during the 
first coupling pulse, which creates significant multimode dy- 
namics, as can be seen in (b). (c) shows that even with large 
amounts of classical noise, it is possible to number squeeze be- 
low the quantum limit. Parameters: (a),(c) and (e): Ati = 1 
/is, At2 = 0.5 /is, thoid = 16 ms. (b), (d) and (f): Ati = 6 
/is, At2 = 0.5 /is, thoid = 3.8 ms. In both cases we assumed 
llo = 50 rad s~^, A^o = 10^, and a 500 Hz spherical harmonic 
trap. 

that when atoms are created in state |2), they are no 
longer in a motional eigenstate of the system. These dy- 
namics are relatively insignificant in first case, when only 
0.25% of the atoms are transferred in the first coupling 
pulse. However, in the second case, ~ 9% of the atoms 
are transferred in the first coupling pulse, and the per- 
turbation to the dynamics during the hold time is signif- 
icant, even though the system is left to evolve for a much 
shorter time. Discussion. — Detection of atoms with 
high quantum efficiency will be required in order to ob- 
serve the squeezing. This is experimentally challenging, 
but has been demonstrated before [3, Furthermore, 
as the parameter space for the experiment is large, a good 
knowledge of the parameters such as trapping frequency 
and Rabi frequency is required such that theoretical mod- 
eling can predict roughly were to look for the squeezing. 
Apart from that, the scheme relies on having good control 
of optical dipole traps and microwave spectroscopy, and 
should be straightforward to implement. An addition ef- 
fect that may degrade the squeezing is atomic loss due to 
inelastic collisions. Using the three-body recombination 



rates recently measured in [l5| , we estimate that roughly 
10% of the atoms from state |1) (the state that we are not 
looking for squeezing) are lost during the 16 ms hold time. 
However, a more important concern is the two-body in- 
elastic collision rate, as it scales as J|?/^o(r)|^ (i^r, that 
is, the same way as Xij ? the nonlinear interaction param- 
eter, so reducing the atomic density will not help, as the 
lifetime due to collisions and the time taken to achieve 
squeezing scale identically. It was observed in [lH that 
with mixtures of different hyperfine states decayed on 
timescales of order several milliseconds. However, if the 
maximal stretched combination of states was used [l6| 
(for example \F = l,mi? = 1), \F = 2,mF = 2)), life- 
times of several seconds were observed. Using this par- 
ticular combination of states would allow for squeezing 
under our scheme, as it has an 7^ ai2. 

It may be possible to create an intensity squeezed atom 
laser via a similar technique to that discussed in this Let- 
ter. By outcoupling two co-propagating hyperfine states, 
and interfering them at particular distance from the con- 
densate, it may be possible to create intensity squeez- 
ing in one of the modes. However, this would require a 
species of atom with two magnetic field insensitive states, 
with scattering lengths an + a22 7^ 2ai2. This require- 
ment could be avoided by using a separated beam path 
interferometer, as the effective Xi2 goes to zero. However, 
it may be difficult to achieve the required mode matching 
in this case. 
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